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Abstract 

This work demonstrates an innovative numerical method for counting and locating eigenvalues with the 
Evans function. Utilizing the geometric phase in the Hopf bundle, the technique calculates the winding of 
the Evans function about a contour in the spectral plane, describing the eigenvalues enclosed by the contour 
for the Hocking-Stewartson pulse of the complex Ginzburg-Landau equation. Locating eigenvalues with the 
geometric phase in the Hopf bundle was proposed by Way [1], and proven by Grudzien, Bridges & Jones [2]. 
Way demonstrated his proposed method for the Hocking-Stewartson pulse, and this manuscript redevelops 
this example as in the proof of the method in [5], modifying his numerical shooting argument, and introduces 
new numerical results concerning the phase transition. 
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1. Introduction 

Way in his PhD thesis [1] , supervised by Bridges at the University of Surrey, developed numerical results 
which supported the hypothesis that parallel transport in the Hopf bundle could locate and measure the 
multiplicity of of eigenvalues for linearizations of reaction-diffusion, differential operators. The eigenvalue 
problem for the operator linearized about the steady state gives rise to a dynamical system on C". The Hopf 
bundle is represented as C C" over the base space CP"“^, and thus has a realization in the phase 

space for an arbitrary dynamical system. By projecting A dependent special solutions in C" onto 
the dynamics induce parallel transport in the Hopf bundle. In the fiber the parallel transport gives a 
winding number that, relative to a reference phase for the contour, counts the multiplicity of eigenvalues 
enclosed by the A parameter, as proven by Grudzien, Bridges & Jones [2]. Way demonstrated his method 
for the Hocking-Stewartson pulse of the complex Ginzburg-Landau equation, utilizing the exterior algebra 
formulation as Afendikov & Bridges did for the Evans function [3]. This work re-examines this example 
and the method— ^ will re-develop this example as in the proof of the method by Grudzien, Bridges & 
Jones [2], and with the example framed in this context, ^introduces new numerical results concerning the 
phase transition. The numerical example highlights two important features that differentiates the method 
of geometric phase from other Evans function formulations; namely, the geometric phase method it utilizes 
either the stable or unstable manifold for the computation of the eigenvalues, and the computation of the 
eigenvalue is continuous in the wave parameter, yielding a continuous accumulation of the eigenvalue driven 
by the system dynamics. 
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2. The Hocking-Stewartson Pulse of the Complex Ginzburg-Landau Equation 


The scaled, complex Ginzburg-Landau equation is given by 

pe^^Yt = Y,, - (1 + iujfY + (1 + iuj){2 + ioj)\Y\^Y 


( 2 . 1 ) 


where p > 0, ip and uj are specified real parameters for the system. The Hocking-Stewartson pulse is the 
steady state solution for the Complex Ginzburg-Landau equation, given by 


Y(x,t) = {cosh{x)) 


( 2 . 2 ) 


Bates & Jones [1] prove that stability of a steady state can demonstrated by stability of the linearization of 
equation (2.1), therefore, consider the linearization about the pulse. 


Definition 2.1. Let C he the linearization of the operator defining equation (2.1) about the Hocking- 
Stewartson pulse pJ 


Considering solutions proportional to e^*, one can derive a non-autonomous system on with asymptotic 
limits in x, as done by Afendikov & Bridges [3]. The system will be of the form 


v' = A{x, \)v n G 

lima._^^oo A(x, A) (A) 

and it is equivalent to the eigenvalue problem in the following sense. 


d 

dx 


(2.3) 


Lemma 2.2. A solution to the A dependent system on is continuous and bounded if and only if it is an 
eigenfunction and A is an associated eigenvalue. 

Proof. This is proven by Alexander, Gardner & Jones and the reader is referred there for the proof. □ 
Remark 2.3. The linearization of the complex Ginzburg-Landau equation about the pulse, C in definition 


2.1: has essential spectrum on the set 

S, 


ess = {p - 1) T 2ip S g 


(2.4) 


and for the parameter values uj = 3, p = ^, and ip = arctan(2) there is a known double eigenvalue at X = 0, 
and simple eigenvalues at approximately A = —6.6357 and A = 15 estimated by Afendikov & Bridges m- 

Definition 2.4. Assume H C C is open, simply connected and contains only discrete eigenvalues of C. 
System p. Jj ) is said to split in LI if A±ao{X) are hyperbolic and each have exactly k eigenvalues of positive 
real part (unstable eigenvalues) andn—k eigenvalues of negative real part (stable eigenvalues), including 
multiplicity, for every A G H. 


Lemma 2.5. For the parameter values uj = 3, p = and ip = arctan(2) system (2.3) splits on the domain 
{A G C : i?e(A) > 0}. Moreover, LI C C can be chosen such that {A G C : Re{X) > 0} C LI and —6.6357 G LI. 

Proof. Afendikov & Bridges [3] demonstrate that the autonomous limits, A±oo(^)) each have exactly 2 stable 
and unstable eigenvalues respectively, for each A such that Re{X) > 0, and in general for A ^ Sess- For a; = 3, 
p = ^, and Ip = arctan(2), the essential spectrum is a curve in C that does not intersect —6.6357; therefore 
an open H C C can be chosen containing A = —6.6357 without intersecting the essential spectrum, and for 
such an Lt, system (2.3) splits on the domain. □ 

Remark 2.6. With the splitting condition satisfied, one may construct the Evans function as done by 
Alexander, Gardner & Jones m — the method of geometric phase was proven by Grudzien, Bridges & Jones 
m under these hypotheses for the Evans function. 
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In order to capture the winding of the unstable manifold of the asymptotic system (A), the compound 

matrix method is utilized, defining a dynamical system on the exterior algebra A^(C^). 


Definition 2.7. Let the matrix A(\, x) define dynamical system of the form (2.3) onC'^. If{zj}1 are vectors 
in €A, the corresponding A^"^) system on A^(C^) = = C® is defined: 


Z' — {X,x)Z X Z 2 ) 

dx 

=A(A, x)zi A Z2 + -zi A A(A, x)z2 
A^^,^(A)= lim A(^)(A,a;) 

ai—>-±oo 


(2.5) 

( 2 . 6 ) 
(2.7) 


Remark 2.8. The system (2.5) yields coordinates for the evolution of the two dimensional subspaces of 
C"^, which allows one to consider the evolution of the unstable manifold for the asymptotic system A_oo(A), 
and particularly the winding it accumulates both in the sense of the relative geometric phase considered by 
Grudzien, Bridges & Jones m and the Chern number of the determinant bundle of the unstable manifold, 
constructed by Alexander, Gardner & Jones m- By way of the proof in m, these winding formulations are 
seen to be equivalent for special solutions which allow trivializations of the determinant bundle. 


Explicitly, Afendikov & Bridges derive the compound matrix system 
Ux = A(A, x)u a: G K A G C 


A(A) 


0 

0 

1 

-1 

0 


032 

0 

0 

0 

0 

0 

042 

0 

0 

0 

0 

1 

—031 

0 

0 

0 

0 

-1 

-041 

0 

0 

0 

0 

0 

0 

-041 

031 

—042 

032 

0/ 


G C® 


( 2 . 8 ) 


with components defined 

031 =Xpcos{'f) + 1 - - (2 - + Sg^) + 6wgig2 

032 = - Apsin(^) - 2 a; - 2(2 - oj)qiq2 + 3 a;(gi + 3 g|) 
041 =Xps\n{'ij}) + 2 uj — 2(2 — u:)qiq2 — — g|) 

042 =Xpcos{ijj) + 1 — uj^ — {2 — uj^){ql + 3 ^ 2 ) — 6 a;gig 2 

and gi, g 2 derived from the expression for the pulse in C^, where 

cos (a; log(cosh(a:))) 


qi = 


g2 =- 


cosh (a;) 

— sin (a; log(cosh(a:))) 
cosh(a;) 


The non-autonomous system (2.8) has the symmetric asymptotic limits 

/ 


Aoo(A) — lima;_^ 4 -oo A(A, x) — 


0 

-p{X) 

7?(A) 

-vW 

-p{X) 

0 


-1 

0 

0 

0 

0 


-p(A) 77 (A) - 77 (A) -p{X) 


0 \ 
0 
1 

-1 

0 

0 / 


(2.9) 
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where the parameters are defined 


p{X) = 2u! + Xpsin{tp) 

77 (A) = 1 — + Xpcos{tp) 


( 2 . 10 ) 

( 2 . 11 ) 


The method of geometric phase makes use of a loop of eigenvectors, chosen analytically in the parameter 
A, that correspond to the dominant unstable eigenvalue for the limiting systems at ± 00 . The method is to 
measure the relative winding the loop at —00 accumulates through its evolution as it meets a corresponding 
loop of eigenvectors at + 00 . An algorithm for constructing these loops of eigenvectors is given by Humpherys, 


Standstede & Zumbrun [0]. However, for the asymptotic system (2.91 such bases ca n be c onstructed explicitly. 
The unique eigenvalues of most positive and most negative real part for system (2.9) are given by (T^,a~ 
respectively, and have associated eigenvectors 


< cr+ = V2^?7 + i/rp +p^, X+(A) = 


M 


< cr = -a/2W?7 + + 


X-(A) = 


2(7+ 

-2p 
(( 7+)2 
-((7+)2 
-2p 

\cr+(((T+)2 - 2ri)) ^ 
( 2(7- 

-2p 

-ia-r 

-2p 

\a-{{a-f - 2p)) ^ 


( 2 . 12 ) 


M 


(2.13) 


These eigenvectors correspond to the Grassmann coordinates for the un/stable subspace of the asymptotic 
system on C^. 

Remark 2.9. For a: < 0 and \ x \ taken sufficiently large, the X dependent initial conditions defined by 


(2.12) approximate the unstable manifold of the fixed point 0 for the asymptotic system, and can be used for 


initial conditions for the numerical method. 


3. The Method of Geometric Phase 

The strategy for determining the eigenvalues of the linear operator C with the Evans function is to choose 
a contour K in the complex plane which does not intersect the spectrum of C and to calculate the winding of 
the Evans function around K. The same strategy applies with the geometric phase formulation calculating 
the winding of the Evans function in the winding in the fibers of the Hopf bundle. This section develops the 
method of geometric phase with C, the linearization of the complex Ginzburg-Landau equation about the 
Hocking-Stewartson pulse. 

Hypothesis 3.1. Assume that the contour K is a smooth, simple closed curve in fl G C such that there 
is no spectrum of C in K. Let K° be the region enclosed by K—assume K° is homeomorphic to the disk 
Z? C and that K is parametrized by A(s) : [0,1] ^ AT with standard orientation. 

The Hopf bundle is realized within the phase space of the exterior system on C^, A^(C‘^) = C®, by 
spherical projection; in general the generic Hopf bundle is realized similarly in C". 

Definition 3.2. The Hopf bundle is a principal fiber bundle with full space 5'^"'-^, base space and 

fiber . The fiber acts naturally on by scalar multiplication by the associated element in the 
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unitary group U(l), and with respect to this action, the quotient is CP" The spaces are related by the 
diagram 

gl -^ S2u-l 


TT 

Cpn-1 

where tt is the quotient map induced by the group action. 


(3.1) 


With respect to this coordinate realization of the Hopf bundle in the phase space C", one may define a 
natural choice of connection. A connection describes the motion of a path within the fibers of the bundle—it 
is possible to define a connection via a connection f-form, as below. For a full discussion on the theory of 
principal fiber bundles and connections the reader is referred to Kobayashi & Nomizu [7]. 

Definition 3.3. For the Hopf bundle viewed in coordinates in C", define the connection 1-form uj 

pointwise for p G as a mapping of the tangent space of the Hopf bundle C Tp{C'^) 


LOp-. rp(5'2"-l) ^ 


Vp ^ {Vp,p)c^ 


(3.2) 


where Q is the Lie algebra of the fiber m- Defining a connection 1-form defines a connection and lo will 
be denoted the natural connection on the Hopf bundle. 

Remark 3.4. The natural connection is a connection of the generic Hopf bundle 5^"“^ and it is the unique 
connection for . This is proven by Way m in §5.5 and the reader is referred there for a full discussion. 


The realization of the Hopf bundle, and its connection, in coordinates for C" also yields an alternative 
formulation for a non-zero path in C", implicitly describing the connection of its spherical projection. 

Lemma 3.5. Let u{t) : [0,1] —>■ C" be a non-zero path, and u(t) be its the spherical projection into 5'^"“^. 
The natural connection of the tangent vector of u{t) is identically 


w — 
\dt 


Proof. By definition 


U) 


d , 


u{t) 


,Im{{u' (t),u{f))) 

* {{u{t)Mt))) 


u{t) 


which verifies equation (3.3). 


(3.3) 


{u{t),u{t)Y {u{t),u{f)y 

u' {t){u{t),u{t)y u{t)Re({u'(t),u{t))) 

{u{t),u{t)) {u{t),u{t))i 

{u'{t),u{t)) _ Re{{u'{t),uit))) 

{u{t),u{t)) {u{t),u{t)) 

, Im{{u'{t),u{t))) 

* {u{.t),u{t)) 


u{t) 


{u{t),u{t))^. 


□ 


With the choice of connection, one may measure the parallel translation of a path in the Hopf bundle—the 
total parallel translation accumulated by a path is referred to as the geometric phase. 
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Definition 3.6. Let 7 C ^ he a path in the Hopf bundle. Then the geometric phase of j with respect 
to the natural connection is the total winding in the fiber given as 




(3.4) 


For a general non-zero path 7 parametrized by u(t) : [0,1] —^ C" one may write the geometric phase as 


GPh) = 


27T 


P Im[{u'(s),u(s))) 
D ((u(s),u(s))) 


(3.5) 


In his thesis [T] , Way developed numerical results which calculated the geometric phase of special solutions 
for the system (2.8). In particular, he considered solutions which corresponded to the stable manifold of 


the system A+oo(A) and, with respect to the equation (3.4), he calculated their phase with respect to the 
contour K. The reformulation of the method of geometric phase, as in the proof by Grudzien, Bridges & 
Jones [2], is elaborated below: 


The Method of Geometric Phase 

• Step 1: Choose a contour K in C that does not intersect the spectrum of the operator C, and let 
A(s) : [0, 1]^ K he a. parametrization of K. 

• Step 2: Varying X G K define a loop of eigenvectors, V=*=(A), for the 44 ^^(A) system 
V^(A) corresponds to the dominant eigenvalue of positive real part. 

• Step 3: Choose xq “close” to —00 and let X~(X,Xo) be a loop of initial conditions. 

• Step 4: Integrate the loop of eigenvectors, following Y' = (A, x)Y, to xi “close” to +00 

• Step 5: Measure the relative geometric phase of X“(A,^i) and V+(A), ie: 

GP(x-{K,xi)^-GP{x+{K)^ (3.6) 

where GP{^) is defined as in equation (|3.5[). 


2.9) where 


Theorem 3.7. Let V"''(A,t) be a solution to the system (2.8) which is in the unstable manifold for the 
asymptotic system A_(„(A). As xi —>■ 00, the relative geometric phase given by (3.6) converges to the 
multiplicity of the eigenvalues enclosed by the contour K. 

Proof. This theorem is proven in a general formulation by Grudzien, Bridges & Jones [2] and the reader is 
referred there for a discussion of the general method. □ 

Remark 3.8. In the above case where T_oo(A) = Ago (A) = A+oo(A), one may take X~{\) = X+(A) so 
that the formulation in equation (3.6) is the difference of the geometric phase of the evolved solution and its 


initial condition—numerically this is given by 


GP{X-{X,xi)) - GP{X-{X,xo)) 


(3.7) 


Note that the relative geometric phase need not be an integer, and indeed cannot always be an integer 
value, as can be seen directly in the proof of the method of geometric phase in [2]. The phase calculation 
must be continuous in a; S (— 00 , 00 ) by the properties of the flow. Moreover, when V“(A) = V+(A) the 
relative geometric phase will always start at zero and transition to the multiplicity of the eigenvalue, as 
demonstrated in the following examples. 
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4. Numerical Results 


This section describes the implementation of the numerical method and demonstrates the transition of the 
geometric phase. In many Evans function formulations, solutions corresponding to the un/stable manifolds 
for the systems Azpoo(A) are integrated to some matching value, usually a; = 0 where the winding of the 
Evans function is calculated [8]. When computing the winding of the Evans function with the geometric 
phase in the Hopf bundle, one considers the unstable or stable manifold and computes the geometric phase 
accumulated relative to the opposite asymptotic condition. The accumulation of the geometric phase thus 
introduces a new dependence on x\ where xi is the final x value of the forward integration described in 
the method of the geometric phase above. As in the examples shown by Grudzien, Bridges & Jones [2], 
the geometric phase undergoes a transition as xi grows large, but for the computation of eigenvalues of 
the linearization of the complex Ginzburg-Landau equation about the Hocking-Stewartson pulse below, the 
transition isn’t uniform across the eigenvalues or the integration parameter. 

In each example below, the contour K is chosen to be the circle of radius .1 about Aq where Aq € 
{0,15, —6.6537}. The contour is discretized into 10,000 even steps, and for each fixed A in the discretization 
of K, the the unstable eigenvector (2.12) is integrated from Xq = —10 forward to some xi. The Matlab 


ODE45 solver is used to find the trajectory of the initial condition A+(A,a;o) with respect to the system 


(2.8), and the trajectory is stored at step sizes of .04 in x. To compute the relative phase in equation (3.6), 
the derivative of a path -^V{X{s)) is approximated via the difference equation 

V {X{s) + Ss) - V {X{s) - 6s) 

Ws ■ 


From the connection equation (3.3), the geometric phase of A+(A) and A {X,xi) is computed with the Euler 


method. The difference of phases in equation (3.6) is computed for each stored value of x and the plot of the 
accumulated geometric phase A“(A,xi) is given for each of the three contours below—because the system 
is symmetric, the relative geometric phase is describe by subtracting the initial geometric phase from the 
terminal geometric phase. However, in each of these examples the initial geometric phase is approximately 
zero, so that the terminal value in the phase plot is approximately eigenvalue multiplicity once the transition 
has terminated. 

The first figure demonstrates the phase transition for the simple eigenvalue at Aq = 15. This figure 
describes a relatively uniform phase transition, exhibited for other examples demonstrated by Grudzien, 
Bridges & Jones [2]. 
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Figure 1: The phase transition plotted for the simple eigenvalue at A Ri 15 


However, the other two plots for Aq € {0, —6.6357} deraonstrate a non-uniform transition both in terms 
of the fluctuation in the phase calculation, as well as the value of x for which the transition begins. 







Geometric Phase 



Figure 2: The phase transition plotted for the double eigenvalue at A = 0. 
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Figure 3: The phase transition plotted for the simple eigenvalue at A —6.6357 

The scale in x direction for the plot of the phase transition at Aq = —6.6357 is longer, ending at Xi = 14. 
Noticeably, the transition here begins later, and doesn’t terminate until it is nearly at the end of the other 
plots, at xi « 10. This example in particular highlights the importance of understanding the phase transition 
for applications. 

5. Discussion 

Evans function calculations are often useful as a stability index [2, describing the multiplicity of eigen¬ 
values of positive real part by computing the winding of the Evans function along the imaginary axis, and 
bounding the integral of the winding along a semi-circle of radius r, as r —>■ oo. In particular, in order 
to utilize the method of calculating the winding with the geometric phase in the Hopf bundle, it will be 
critical to understand the nature of the phase transition. As demonstrated in the example above, the phase 
transition is neither uniform in the integration of the x direction, nor uniform across eigenvalues—indeed the 
calculation may fluctuate and the initiation and termination of the transition differs for each of the above 
results. Eor utilization as a stability index, one must understand the relationship between the transition and 
the underlying steady state to efficiently compute the eigenvalues. The geometric phase must eventually 
converge to the multiplicity of the eigenvalues enclosed by the contour, but a theoretical understanding of 
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the transition of the phase will be an important development for both the numerical method and the un¬ 
derstanding the eigenvalue problem itself—indeed the method of the geometric phase offers a unique insight 
into the continuous accumulation of the eigenvalue as driven by the system dynamics, a new insight not 
afforded by other Evans function methods. 

Currently the method of geometric phase is limited by the dependence on the exterior algebra formulation— 
for usual systems on C", where the stable and unstable manifolds are of dimension approximately the 
dimension of phase space for the exterior algebra grows approximately exponentially in n, as discussed by 
Humpherys & Zumbrun [9]. However, the fact that the method of geometric phase relies only on either 
the unstable or stable manifold for the eigenvalue calculation highlights the potential for future reductions. 
Reducing the calculation of the phase to a frame of solutions spanning the un/stable manifold, rather than 
the wedge product, is the subject of active research. 

6. Conclusion 

Way developed numerical results for the calculation of the geometric phase for solutions in the stable 
manifold of the H_^(A) system for the linearization of the complex Ginzburg-Landau equation about the 
Hocking-Stewartson pulse, in his thesis [T]. This work builds on that discussion, framing the numerical 
method in terms of the relative phase as in the proof of the method of geometric phase by Grudzien, Bridges 
& Jones [2], as well as demonstrating the phase transition in the calculation with new numerical results. In 
studying the numerical method, understanding the phase transition and its relationship to the underlying 
wave will be critically important for finite approximations and particularly for use of the geometric phase for 
stability indices. The varied nature of the phase transition across eigenvalues for a single example highlights 
the need to understand this transition, and this work opens new theoretical questions for the advancement 
of Evans function techniques, and the method of geometric phase. 
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